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Abstract 

Motivated by anomalies in cosmic ray observations and by attempts to solve questions of the Standad Model of particle 
physics like the (g — 2) M discrepancy, (7(1) extensions of the Standard Model have been proposed in recent years. Such U(l) 
extensions allow for the interaction of Dark Matter by exchange of a photon-like massive force carrier 7' not included in the 
Standard Model. In order to search for 7' bosons various experimental programs have been started. One approach is the 
dedicated search at fixed-target experiments at modest energies as performed at MAMI or at the Jefferson Lab. In these 
experiments the process e(A, Z) — > e(A, Z)l + l~ is investigated and a search for a very narrow resonance in the invariant mass 
distribution of the l + l~ pair is performed. In this work we analyze this process in terms of signal and background in order to 
describe existing data obtained by the Al experiment at MAMI with the aim to give accurate predictions for exclusion limits 
in the 7' parameter space. We present a detailed theoretical analysis of the cross sections entering in the description of such 
processes. 
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I. INTRODUCTION 



With the recent observation of a new boson at the LHC, which is expected to be the Higgs boson, the last missing 
element of the Standard Model of particle physics (SM) seems to be discovered 0, 0|- Despite of this, nowadays the 
existence of dark matter, which is not included in the SM, is established as a necessary ingredient in order to explain 
the energy density of the universe within the cosmological standard model ■ The nature of dark matter is however 
still a wide open question. Neither is it known what dark matter is made of, nor in which way it is interacting with 
other particles, e.g. the SM particles. Besides the unsolved problem of dark matter, the SM itself does contain several 
issues like the discrepancy in the theoretical and experimental determination of the anomalous magnetic moment of 
the muon (g — 2) M , the proton radius puzzle, or weak scale questions like the hierarchy problem, which all could be 
hints for physics beyond the SM. 

Recent observations of anomalies in astrophysical data [fUll have motivated to consider extensions of the SM by 
including an additional U(l) gauge group which could explain such anomalies (^, Though the idea to extend the 
SM by an additional U(l) recently became popular, it did not rise up with the observations. In many well motivated 
SM extensions, e.g. from string theory, additional U(l) groups appear naturally (TU - H^ . 

Extending the SM by such an U(1)d group generates an additional gauge boson 7' which is able to interact with the 
electromagnetic current of the Standard Model. Although this interaction is forbidden at tree level it is possible via 
kinetic mixing [l3| giving rise to an effective interaction Lagrangian 

Ant = i ee ^ S M 7 M 4>sm A'^ , 

where A' denotes the 7' field. Furthermore, e is the kinetic mixing factor parameterizing the coupling strength relative 
to the electric charge e, and describes the interaction of the additional gauge boson with the electromagnetic current. 
The 7' may gain a mass m.y which can be estimated to be in the range of lOMeV to a few GeV 
mixing factor e 2 = a! /a is predicted from various models to be in the range 1CP 12 < e < 10~ 2 
coupling via kinetic mixing the 7' may decay to dark matter particles as well as SM matter particles. In the case, 
that the decay to dark matter is kinematically forbidden and my > 2m e — 1.022 MeV, which this work will focus on, 
the 7' will decay to SM particles and therefore must be observable at accelerator experiments. 

The 7' interacts with SM particles and has properties which are very similar to that of the photon. Since by now 
such a boson could not be observed one often refers to the 7' (which is also denoted as A' , U, 4>) as heavy, hidden, 
para- or dark photon. Within this minimal model the free parameters are the mass my and the coupling strength 
e. In pioneering works several constraints from existing data were obtained on these parameters e.g. in beam dump 
searches or by the BaBar experiment, as well as from (g — 2) analyses (TBI. l2lj. 

The coupling of the 7' to SM particles and the predicted mass range allows for the 7' search by accelerator experiments 
at modest energies with high intensities. While collider experiments are ideally suited for higher 7' masses, fixed-target 
experiments with their high luminosities are ideally suited for the 7' search in the MeV to 1 GeV range [l9l. [2l| - [2^ | . The 
proposal to search for the hidden gauge boson by fixed-target experiments motivated several experimental programs, 
both by the Al collaboration at the MAMI accelerator in Mainz [26| as well as at the CEBAF facility at Jefferson Lab 
with the APEX (25l[27j. HPS [28| and DarkLight [^HH experiments. The Al and APEX experiments already have 
published first data. Furthermore in many recent publications, constraints on the 7' parameter space from the analysis 
of beam dump searches [3ll - l34j | , meson decays and collider experiments (35l - [38| as well as from other arguments were 
given [3^. l4(il|. and are summarized in Fig. 1121 In addition, many other experiments were proposed to probe the light 
hidden sector or are underway, for a review see e.g. Ref. [4l|. 

In all considered fixed-target experiments an electron beam is scattered off a fixed target which is either a proton 
or a heavy nucleus like tantalum. Induced by this electromagnetic process a 7' may be radiated from the electron 
beam and decays into SM particles like an electron-positron pair. Detecting the decay particles and reconstructing 
the invariant mass of the pair allows to search for the hidden gauge boson by a bump hunt. The 7' will manifest itself 
by a very sharp peak above the radiative background that results from the corresponding process where a virtual 
photon is radiated from the electron beam which creates a lepton pair, too, i.e. the underlying process 

e(A, Z) e(A,Z)l + r 

is investigated. 

If there is no bump seen in the invariant mass spectrum this allows to exclude regions of the 7' parameter space given 
by the kinetic mixing factor e and its mass my . In order to perform this study a precise knowledge of the signal and 
background cross sections are crucial. Such precise study is the main subject of the present work. 
This work is structured as follows: In section [TT] we present our calculations of the signal and background cross 
sections. In section IIIII we present our results of the cross section calculations for the experiments performed at 
MAMI. Furthermore we present a comparison with available data. In section IIVI we propose new searches at MAMI 
and at the new MESA accelerator and present our predictions for the exclusion limits. 
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Figure 1. Tree level Feynman diagrams contributing to the ep — > epl + l~ amplitude. Upper panel: exchange of the timelike 
boson V and a spacelike 7 (TL). Lower panel: the spacelike boson V and a spacelike 7 (SL). In addition to these direct (D) 
diagrams the exchange term (X), which consists of the same set of diagrams with scattered electron and electron of the e + e~ 
pair exchanged, also contributes. 



II. CALCULATION OF THE SIGNAL AND BACKGROUND CROSS SECTIONS 

The underlying diagrams for all fixed target experiments mentioned so far are shown in Fig. [TJ We calculate this 
process exactly in leading order of QED and furthermore apply leading order radiative corrections of the corresponding 
elastic scattering process to obtain an estimate of these corrections. 

An electron beam of energy Eq is scattered off a fixed target, which may either be a nucleon or a heavy nucleus of 
atomic numbers (A, Z). In the following, the target mass M refers to the nucleon mass Mn or to the mass of the 
heavy nucleus Ma -ix Mn- As subprocess to the elastic scattering an intermediate vector particle V is produced 
and creates a lepton pair (l + where the lepton mass is denoted by mi. Although the existing and planned fixed 
target experiments only consider electron-positron pairs in which a bump hunt is performed, our calculations are 
performed generally for any kind of lepton species (i.e also applies to the case), i.e. we do not neglect the mass 

of the lepton. 

The isolated 7' production process is given by the coherent sum of diagrams (a) and (b) while the background, 
resulting from the exchange of a virtual photon, is given by the sum over all diagrams, where the intermediate vector 
particle V in diagrams (a) and (b) is 7' and 7*, respectively. 

We assign a finite decay width Ty to the 7'. The partial decay width to a SM lepton pair is given by 

T ^ l+l - = ^,yJ m Y~ 4m i K' + 2m f)> 

with a = e 2 /(4n) ~ 1/137. 

For kinematically forbidden decays to dark matter the total width can be related to the partial width by Fy = 
JV e ff x r 7 /_j.;+;-, where iV ff is a weight to account for other degrees of freedom in SM decays. Since in this case 
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the width is very small, only a small mass window around the peak will contribute to the signal and thus the cross 
sections for real and virtual 7' multiplied by N e g are equal (2lL l29l|. 

In this work we denote the four-momenta of the initial and final beam electrons by k — (Eq, k) and k' = (E' e , k')\ the 
four-momenta of the initial and final target state by p = (E p , p) and p' = (E' p , p') and the lepton pair four-momenta 

by Z_ = [E-, 1-) and l + — (E+, Z + ), for the lepton and anti- lepton, respectively. The initial and final electron spins 
are denoted by Sk and s' k ; the spins of the initial and final proton by s p and s' p ; and the spins of the created lepton 
and anti- lepton by s_ and s+. Furthermore we follow the conventions of Bjorken and Drell [42l |. 
The invariant amplitudes required to calculate the cross section can be read of from these Feynman diagrams. As 
in the two diagrams in the upper panel of Fig. [T]the intermediate boson V is timelike, we refer to this amplitude as 
TL. Correspondingly, we refer to the diagrams in the lower panel, where the V is spacelike, as SL and their sum is 
denoted by SL + TL. 

In the case that the l + l~ pair and the beam lepton are of the same species as for the existing experiments, another 
set of diagrams is allowed. Since one cannot distinguish the electrons in the final state, the same diagrams of Fig. [T] 
with the scattered (beam) electron and created electron of the pair exchanged, also have to be taken into account. 
Therefore, following the notation of Ref. [42j, we refer to the diagrams depicted in Fig. [T] as "direct" contribution and 
to those with exchanged final state electrons as "exchange" contribution, labeled by D and X, respectively. 
For the TL diagrams one finds for the isolated 7' production process 
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where Ty denotes the total 7' decay width. The amplitude of the 7* background is given by: 

7 p 4 —n a P 

LfTL _ l _± y j -pair icy, 

[P ~ P) 1 

where the external momenta are denoted by q = k — k' , q' = Z_|_ + I— as on Fig. [1] Furthermore the leptonic tensors 
are given by 

3-m = u e (k',s' k ) ( 7 P ^ ■ J - + 9 m 7a + j a ^ j k + J 2 - + ) u e {k,s k ), 



(k-q'Y-m 2 (k' + q'Y 



m 



with m denoting the mass of the electron. While in the case of a proton target the hadronic current J^, is given by 

= un(p', s' p ) r M u N {p, s p ), 

with the parametrization of T fJi (Q 2 ) = Fi(Q 2 ) 7^ + F2(Q 2 ) i a^q" flM using the Dirac and Pauli form factors F\ and 
F2 and Q t = — (p — p') 2 > 0. For a heavy nucleus it can be written to good approximation as 



J»=Z-F(Q t )-(p + p' 



where F(Q t ) = 3/(QtR) 3 ■ (sin (Q t R) — QtR cos (Qt R)) is the nuclear charge form factor with R = 1.21 fm • Ai . 
The nucleus spin as well as contributions from the breakup channel and nuclear excitations can be neglected to good 
approximation. Effects due to the nucleus spin are suppressed by the large nucleus mass, which can be checked 
analytically. The inelastic contribution can be neglected since the momenta transfered to the nucleus are small. 
The numerator of the 7' propagator in Eq. ([TJ) can be simplified as (—g a ^) since the four-momentum q' is contracted 
with the lepton current jP air and thus the second term vanishes due to current conservation. 
For the SL diagrams the invariant amplitude is given by 

IP -P) 1 



with 



,| (7 • (g - l+))+mi (7 ■ (j- - q)) + mi , 

X m = Ul(l-,S-) 7 M — —2 — 7a + 7a— "2 T~ 7m "K«+. S +)i 

[q — l + j — mf («_ — q) — 1 

j h ^ = u e (k',s' k )^u e (k,s k 
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Figure 2. Tree level Feynman diagrams of the double VCS contribution. 



Although the virtual 7' exchange via the SL process is not forbidden, it will not be considered here as it would not 
result in any bump in the e + e~ mass spectrum. The propagator in Eq. ^ in that case would be replaced by 



-g aP -g af) 



9 9 2 ' 

q- q A - m*, 

and due to the spacelike q 2 < for scattering processes the denominator always leads to a suppression of this 
contribution, whereas the denominator in Eq. ([T]) leads to a peak in the signal. Thus this contribution of virtual 7' 
exchange via the SL process to the cross section can be neglected. 

In the case of a proton target another important contribution, the double virtual Compton scattering (VCS), 
emerging from the third set of Feynman diagrams shown in Fig. [3J appears. In case of a heavy nucleus target this 
term is strongly suppressed due to the large mass. In this work we will restrict our study to estimate the influence 
of the nucleon pole contribution drawn in Fig. [5] which serves as a good approximation. The invariant amplitude is 
given by 

v CS _ -ie*e 2 -9 p + q q p /m y lr 

t' ~ q 2 ^-m^+zmyr/^^^ ' () 
for the isolated 7' production process, and 



M VCS _ 16 y -li nj -pair /j-x 
JVl 1* — 2 gl2 Acam J p \p) 



for the 7* background, with 



n, - u M r t , n 7- (p-q' ))+MN r , , r , a 7- (g ' + Q ')) + M N ^ ,\ 

7~tfj,a =u p (p,s ) T^qt+q)— ——T a (-q) + T a (-q)—— -5 ——T^qt+q) u p (p,s p ). 

V ip-q') - Mn {p' + q') - m 2 n J 

As mentioned before, the electron from the scattered beam and the one from the lepton pair cannot be distinguished 
and besides the direct term the exchange term has to be accounted for. Therefore the full amplitude of the process 
reads as 

My+r = (My L + + My) - {(M^ + + M^) (e~ o V )) 

= (Ml L Y + Ml h r + Mg>) - (Ml L Y + Ml L r + X| L 7 .) , (6) 

for a heavy nucleus target and 

M Y+r = {Ml L y + Ml L r + Ml L r + M^%) - (Ml L Y + Ml L r + X| L 7 , + M%%) , 

for a proton target. In the second term of Eq. ([6]) all quantities associated with the scattered electron and the pair 
electron are exchanged. The exchange 7' term can be neglected, as the 7' propagator does not peak and thus a 
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possible signal is suppressed by e 2 . Due to the exchange of final state electron momenta, the amplitude describing 



the signal A4^, L , as well as the background contributions A4]y h 



7* 



and 7W| L 7 , 



contain a structure 



(7 • (fc - I- - l+))+m 
(k-1-.-l+f -m? 1 

that contributes to the irreducible background. This leads to a large contribution from A^x L 7 * m the case of forward 
scattering, since the denominator of the propagator is close to zero. Forward scattering was proposed to enhance the 
signal strength, while not increasing the background A4zk. Taking the background contribution A4y^ » into account, 
this argument is not applicable anymore, since now the background is also enhanced. 

The cross section of the ep — > epe + e~ process is computed from the general expression for 2 — > 4 particle processes 



da = 
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4-y/ (k ■ p) 2 — m 2 M 2 
d 3 k' d 3 p' 



(2tt) 4 <5 (4) (k + p-k' -p' -1-- 1+) 

d 3 C d 3 C — 
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Using a convenient set of variables we can express the cross sections as 
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(8) 



pair, 



= (p + q) is the Mandelstam invariant of the 7*- 
— (M + mii) 2 ^ — (M — m;/) 2 ^ denotes the kinematical triangle function. 



where mu = y l' 2 1S the invariant mass of the I 

target subprocess, A(s, M 2 , mf^ 

Furthermore, we use the superscripts Q L to label the lab frame, ()* for the (q + p) rest frame, and ()** for the q' 
rest frame. This approach avoids ambiguities in the kinematics and automatically gives the full kinematically allowed 
region of the phase space. In the considered type of experiments only a small fraction of the kinematically allowed 
phase space is probed. The allowed region is given by the detector acceptances in the lab frame. Therefore it is 
convenient to calculate the cross section directly in terms of lab frame quantities and to use the recursively built up 
phase space as a cross-check. 

Since fixed target experiments are considered here, the target four- momentum p simplifies to p — (ill, 0). Furthermore 
in the considered experiments the detectors and the beam are aligned in the same plane which we account for by the 
choice of our parametrization of the momentum vectors of the detected particles. Since neither the scattered hadron 
nor the scattered electron will be detected in the experiments, as long as the electrons are treated as distinguishable 
particles, the dependence of the cross section on their four-momenta has to be eliminated. Therefore the three- 
momentum conserving (^-function is used to eliminate the three-momentum of the final hadron state p' and energy 
conservation is used to express the absolute value of the three-momentum of the scattered electron \k'\. The remaining 
dependence of the cross section on the electron scattering angle is removed by integration over the full solid angle fi e ' . 
Furthermore one is interested in the cross section as function of the invariant mass of the created lepton pair, which 
is equal to the squared four- momentum of the intermediate vector boson q' 2 = m 2 ; . Therefore we trade the absolute 
value of Z_ for q' 2 . 

Thus one finds from Eq. ([7]) for the differential cross section in the lab frame 
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d8 t 
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(9) 



where this equation is understood to be evaluated with \l + \ and \k'\ given in Eqs. (|A5[1 and (IA6[) . and 



9<5i 
d\k'\ 



and 



are given by Eqs. (|A7|1 and (|A8[) . respectively. A more detailed derivation of the cross section is presented in 
Appendix fA"! 

Furthermore, we will apply radiative corrections of elastic electron-proton scattering to the cross section to achieve a 
better comparability with the experimental data. Therefore the cross section of Eq. © is multiplied by Eq. (A71) of 
Ref. . By applying these radiative corrections the value of the cross section is reduced by an amount in the range 
of 10 - 20 %. 

The comparison with experimental data can be performed by integrating Eq. ([9]) over the experimental acceptances. To 
obtain the acceptance integrated cross section Act, which can be related to experimental count rates by multiplication 
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with the luminosity, a non-trivial 8-fold integration is necessary. Furthermore the structure of the squared matrix 
element contains several strongly peaked structures which makes the numerical calculation of this integral challenging. 
Any of the fermion and photon propagators in the Feynman diagrams shown in Fig. Q] can possibly be near the mass 
shell in a certain kinematical setting. Although there is no real divergence existing, since the non- vanishing mass of the 
electron serves as a regulator, the calculation of this strongly peaked structures either needs further approximations 
or a large numerical effort. In our study we try to use as less approximations as possible. We thus decide to use an 
integration method that allows to deal with these peaked structures by increasing the numerical precision. Therefore 
for the numerical integration the VEGAS algorithm [44J has been chosen, which is a well established Monte Carlo 
integration method in particle physics. The standard deviation and the x 2 of the result of the integration are used to 
decide whether the computed value is reasonable or not. During our calculations it turned out, that - at least for the 
case of MAMI kinematics - one cannot use a vanishing electron mass to achieve numerically stable results. 
In order to perform these calculations in a reasonable amount of time, we have performed a highly parallelized 
calculation. Therefore the integral is computed on Graphics Processing Units (GPUs) using the NVIDIA CUDA 
framework [45| and the implementation of the VEGAS algorithms on GPUs published in Ref. j46j. The use of the 
GPU version reduces the time needed for the evaluation of the acceptance integrated cross section by a factor of ~ 60. 
We have checked the results achieved by the GPU calculation with ordinary calculations on CPUs and find that for a 
same numerical precision the results are equal within their standard deviations, which are below 10~ 4 relative to the 
obtained value. 

The radiative background is described by the acceptance integrated cross section 

Aa 7 cx \(M^+M^) - ({M^+M^) (e- O , (10) 

where the prefactors on the right-hand side are the same as appearing in Eq. ©. For later use, besides the cross 
section of the process including 7* and 7', we define the direct timelike radiative background cross section and the 
direct timelike 7' cross section as 

+ AfeA 7 .| 2 , (11) 
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(12) 
(13) 
respectively. 

In order to compute exclusion limits on the coupling strength e from existing data, a relation between the cross 
sections of Eqs. (TTU1) and (fT3")) giving rise to e is required. We split the 7' + 7 cross section as 

Aoy +7 = Aer 7 + Aoy + Aer int , 

with A<7int denoting the interference part. Dividing Eq. ([TTj) by Eq. (TTU1) leads to 

Aoy +7 ^ 3ir e 2 my A<r 7 L Aa- ln t 
Aer 7 2N a 8m Aer 7 Acr 7 

We have used Eq. (19) of Ref. [21| in order to approximate the ratio of ay and ct 7 l as 



Aay 3ir e 2 my 



where N is the ratio of the decay widths r 7 /_j. e + e - and r 7 /_>. M + M - taking other possible final states into account and 
Sm is the experimental mass resolution, i.e. the mass bin width. For 7' masses > 400 MeV hadrons also contribute to 
the final state and thus our parametrization of N is not valid anymore. In the 7' mass range considered in this work 
only electrons and muons are contributing as possible final states. Our numerical calculations for a wide range of 
parameters my and e of the interference part a- m t from the cross sections (fTU)) . (fTTj) . and (fT3")) show that the interference 
between 7' signal and QED background can be neglected. We find, that A<7i n t/A(7 7 is less than 10 -3 , which is in the 
range of the achieved numerical precision. Furthermore we find a very good agreement of the approximated ay /<r 7 L 
with our exact calculation for the largest part of the parameter region for my and e. Therefore e can be computed 
from the cross section ratio as 

a = / Agy ±7 \ A<7 7 2Na Sm 

\ Aa 7 J ActT l 3^ m Y ' K ' 
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Table I. Acceptances of the used spectrometers A and B at MAMI (47T |. 
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Figure 3. Comparison of theory calculations and experimental data for a m e + e - bin width of 0.125 MeV. Black points: Data 
taken in a particular run of the MAMI 2010 experiment 26] in setup 1. Solid curve: Theory calculation of the background 
cross section. Dotted curve: Theory calculation of the background cross section without radiative corrections. Dashed-dotted 
curve: Theory calculation of the direct SL + TL cross section. Dashed curve: Theory calculation of the direct TL cross section. 

The ratio Aer 7 ' +7 /A<7 7 is the (aimed) signal sensitivity, which has to be determined from the experiment. Furthermore, 
by using the ratio Acr 7 '+ 7 /Acr 7 for the extraction of e 2 , possible effects not accounted for in our approximation of the 
nuclear current will cancel each other. For the prediction of exclusion limits we estimate AtJ 7 '+ 7 /AtT 7 — 1 as signal 
over background ratio 




#5 ^/Aa y x V 



where and are the numbers of signal and background events in one mass bin, respectively, and L is the 
integrated luminosity. The factor 2 results from the fact, that in agreement with other publications we determine the 
exclusion limits on the 2 a level. Since the exclusion limit on the coupling strength e 2 is depending linearly on the 
ratio of the background cross section Acr 7 to the TL cross section with distinguishable final state electrons Act 7 L , the 
precise knowledge of these quantities is crucial to obtain an accurate result. Therefore the next section will deal with 
the analysis of these background ratios for the existing experiments. 

III. COMPARISON OF EXPERIMENTAL DATA AND THEORY CALCULATIONS FOR MAMI 

Two dedicated fixed target experiments, one by the Al collaboration at MAMI (2|| and the APEX experiment at 
JLAB [27], have already started taking data. 

A. Test run 2010 

A first test run to proof the feasibility of a dedicated 7' fixed target search experiment has been performed at MAMI 
by the Al Collaboration in 2010 [26]. In this experiment no evidence for the existence of the 7' could be found and an 
exclusion limit on the 7' parameter space was formulated. A sample of the data taken in this experiment compared 
to our calculations can be seen in FigJ3] 

The kinematical settings of this experiment can be taken from Table 1 in Ref. [26j. For the comparison of the calcu- 
lation and the data, the setup 1 as given in Ref. [26] was chosen, since for this setup a luminosity measurement has 
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Figure 4. Angular distribution per 0.5° with respect to the polar angle of the scattered electron for the MAMI 2010 experiment. 
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Figure 5. Calculated direct (left panel) and exchange (right panel) term of the cross section assuming distinguishable electrons 
in the final state. Solid curve: SL + TL cross section. Dashed curve: TL. Dashed-dotted curve: SL 
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Figure 6. Solid (dashed) curve: Ratio of the background cross section Ao" 7i d+x (Act 7i d) to the direct TL cross section Acr-J 
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Table II. Kinematics of the MAMI 2012 7' search. Electron scattering angle: <j>- — 20.01° (spectrometer A). Positron scattering 
angle: 4>+ = —15.63° (spectrometer B). The number in the label of the kinematics refers to the invariant mass around which a 
setting is centered. 

been performed, finding an integrated luminosity of C = 41.4 fb -1 for the selected sample of events. A background 
contribution of around 5% was already subtracted in this sample, the systematic uncertainty in the luminosity from 
the knowledge of the thickness of the target foil is below 5%. The acceptances as shown in Table U have been used 
as integration limits for the theory calculation. Unless mentioned otherwise, the mu integration is performed over a 
range of 0.5 MeV, which is equal to the typical FWHM mass resolution of the considered experiments. 
As seen on Fig. [3J our calculation (solid curve) of the radiative background given by Eq. (|10p and the experimental 
data (points) are in good agreement. Due to our estimate of the nuclear current and of the radiative corrections we 
expect the small discrepancy between theory and data seen from Fig. [3j The influence of the radiative corrections is 
displayed by the solid and dotted curve on Fig. [3] which are calculated with and without radiative corrections, respec- 
tively. It is obvious from Fig. [3J that the applied radiative corrections lower the result of the theory calculation by 
an amount in the range of 10 — 20 %, as mentioned in section[II] The calculation of the full QED radiative corrections 
for such a process is very involved. However, one can see from Fig. [3J that our approximate treatment of the radiative 
corrections already provides a very good approximation, as theory and data already are in good agreement. 
The dashed (dashed-dotted) curve shows the direct TL (SL + TL) cross section. This indicates, that a large contri- 
bution to the cross section results from the antisymmetrization due to the indistinguishability of the scattered beam 
electron and the pair electron. The kinematical setting has been optimized to reduce the SL background. 
The angular distribution with respect to the polar angle of the scattered electron presented in Fig. 2] points out, that 
for the 2010 Al experiment the crossed TL amplitude is responsible for a second peak in the background cross section 
compared to the direct amplitude (dashed curve) which only peaks at very forward scattering followed by a rapidly 
dropping tail. The exchange SL term nevertheless enhances the tail of the angular distribution significantly. 
Fig. [3] reveals, that in the chosen kinematic setting the exchange term contribution is about twice as large as the direct 
SL part, which initially should be minimized. This means, that the largest contribution to the radiative background 
does not originate as assumed from the processes given by the direct SL Feynman diagrams of Fig. [1] but from the 
processes described by diagrams with exchanged final state electrons. 

For the investigated kinematic setting we calculate the ratio of the background cross section to the direct TL cross 
section which is the crucial quantity entering the determination of the exclusion limit on e 2 , according to Eq. (|14j) . 
One notices from Fig. [6] (solid curve) that the ratio Act 7! d+x/Act^ l smoothly varies between 15 and 25 for most of 
the invariant mass range. Neglecting the necessary contribution of the exchange term to the cross section, the ratio 
is lower by a factor of about 3 for the investigated range (dashed curve on Fig. [3J). 



B. 2012 

The Al Collaboration started a 7' search run at MAMI in 2012, probing the kinematics given in Table ILT1 in which 
no signal of a 7' was found. The obtained invariant mass distributions can be seen in Fig. The invariant mass 
distributions calculated from the different cross sections are compared: background (solid curve), SL + TL exchange 
term (dotted), SL exchange term (double-dashed), SL + TL direct term (dashed), and TL direct term (dashed-dotted). 
It turns out that the SL exchange process is the largest contribution to the radiative background. Fig. [7J illustrates 
the dependence of the seperated background contributions on the invariant mass m e + e -. At low invariant masses the 
SL exchange term dominates the cross section. Although the SL direct and the TL exchange terms become more 
important for increasing m e + e - , the SL exchange term remains the largest contribution to the cross section. The ratio 
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Figure 7. Simulation of the invariant mass distributions calculated from the different cross sections for the kinematics probed 
at MAMI in 2012: background (solid curve), SL + TL exchange term (dotted), SL exchange term (double-dashed), SL + TL 
direct term (dashed), and TL direct term (dashed-dotted). 
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Figure 8. Combined plot of our result for the ratios Act 7 /A<t^ l of each setting, starting with the lowest beam energy on the 
left. 
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Figure 9. Invariant mass distributions from the feasibility study for the MESA experiment. Solid curve: SL+TL (direct 
exchange term), dashed curve: direct TL, dashed-dotted curve: direct SL+TL. 



between the TL direct term and the SL exchange term has a similar behavior, retaining nearly the same maximum 
value in each of the considered settings. Furthermore, Fig. [7] shows the importance of the interference parts of the 
cross section, which are necessary to describe the data correctly. 

In Fig. [5] we present a combined plot of our result for the ratio Aa 7 /Acr^ L for each setting given in Table HTl which is 
crucial to obtain the exclusion limits on the 7' mass m 7 < and its coupling strength e 2 following Eq. Q140. as function 
of the invariant mass m e + e - . Due to the particular choice of kinematics in that experiment, the ratio Aer 7 /Ao\^ L has 
a value between 10 — 15 in the probed mass range. 

In Fig. [12] our predictions for the exclusion limits on e 2 for this set of kinematics are indicated by the dashed curve 
for an assumed integrated luminosity of around 10fb _1 . 
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Figure 10. Left panel: Invariant mass distributions from the feasibility study for the MESA experiment for a proton target. 

Solid curve: SL+TL (direct + exchange term), dashed curve: direct TL, dashed-dotted curve: direct SL+TL. 

Right panel: Isolated VCS cross section. Dashed curve: exchange term contribution, dashed-dotted curve: direct contribution. 
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Figure 11. Combined plot of our result for the ratio Ao" 7 /Acr^ L of each setting for the MESA experiment. The settings 
correspond with the following beam energies and scattering angles (from left to right): Eq = 80, 120, 160 MeV with 4>t = ±10°, 
E = 120 MeV with T = ±20° 

IV. FUTURE SEARCHES AND DISCUSSION 

Recently the construction of the Mainz Energy Recovering Accelerator (MESA) has been approved. MESA is aimed 
to provide a high intensity electron beam up to beam energies of about 160 MeV and thus should be ideally suited to 
probe the 7' parameter space for low masses. In this section we perform a feasibility study to carry out this search 
by using two small spectrometers. We assume, that each of these spectrometers has a horizontal and vertical angular 
acceptance of ±50 mrad and a momentum acceptance of ±5%. A possible 7' experiment at MESA can be performed 
using a gas target to minimize the multiple scattering in the target material. Therefore, applying the same program 
code as in section IIII1 we perform our calculations using a Xenon target in order to obtain as large cross sections as 
possible. The integration over the invariant mass m e + e - is performed for a 0.125 MeV interval. 

The results for the obtained invariant mass distributions of this study are shown on Fig. |H1 The kinematics were chosen 
such that the central scattering <p of the e~ (e + ) is ±10° (—10°) and the central momentum is \l \± = 0.98 x Eq/2 for 
beam energies Eq of 20, 40, 80, 120, and 160 MeV. Furthermore, we have calculated one setting for Eq = 120 MeV 
and (f> T = ±20° in order to cover the full so-called (g — 2) M welcome band together with the MAMI 2012 settings. We 
assume a beam time of about 3 months and a luminosity of 10 34 cm _2 s _1 . 

Since the low mass region m 7 ' < 10 MeV in the (g — 2) M discrepancy is already excluded by the electron anomalous 
magnetic moment (g—2) e , the settings for beam energies of 20 and 40 MeV will not enter the exclusion limit calculation. 
Therefore we do not have to deal with the difficulties in the low mass regime. From our exact calculation of the signal 
cross section Aay we find for the considered range of parameters a good agreement with the approximation of the 
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Figure 12. Compilation of existing exclusion limits and our predictions: For a better visualization we restrict ourselves to the 
region currently accessible at fixed-target experiments. Only existing limits as published in Refs. [l6, 21, 36, 37, 39, 4(| are 
shown, represented by the shaded regions. We do not show the predictions for other experiments [25l. I28M3C | in this figure 
which are scheduled to probe the same region of parameter space. The limits of MAMI and APEX are those as given in their 
publications [26l. l27l|. The prediction of this work for the exclusion limit expected for the MAMI 2012 experiment discussed 
in section fill Bl is depicted by the dashed curves. The prediction for MESA obtained in section [TVl is indicated by the dotted 
(dashed-dotted) curve for the setups with a central scattering angle of 10° (20°). 



signal cross section given in Ref. [2lJ. 

For comparison we show in Fig. 1101 the acceptance integrated cross section depending on m e + e - for a proton target 
with a beam energy of Eq = 80 MeV. In the left panel the same curves as in Fig. are plotted. In the right panel of 
Fig-Hniit is demonstrated, that the VCS contribution corresponding with the Feynman diagrams in Fig.[2]are smaller 
by more than 6 orders of magnitude in the chosen kinematic setting, and can thus be neglected. As indicated by the 
shape of the curves for Acr^jt^x and Aa^ u in Figs. 191 and [TU1 the ratio of these two quantities is equal and thus 
the kind of target does not affect the exclusion limit concerning the QED background. 

Fig. rm shows the calculated ratio Acr 7 / Aa^ L which reaches a value around 8 — 10 for the proposed settings. The 
expected exclusion limit on e 2 as obtained from Eq. (fT4")l . to the invariant mass spectra of Fig. [51 is presented on 
Fig.HH where a mass resolution of 0.125 MeV was assumed. The dotted (dashed-dotted) curve on Fig. [T2l represents 
the settings with central angle of 10° (20°). At very low masses below 10 MeV Eq. (|14p does not serve as a good 
approximation for the exclusion limit anymore, since Eq. (19) of Ref. (2l| overestimates the 7' signal cross section by 
up to 50%. 



A compilation of the existing exclusion limits is presented in Fig. 1121 which shows the region 5 MeV < m 7 ' < 
600 MeV and 10 -8 < e 2 < 10 -4 accessible at fixed-target experiments. Furthermore, existing limits as published in 
Refs. (l6l . I21L [H, [I?], OS 4C are also shown, and are represented by the shaded regions. Let us mention that other 
planned experiments |25l . 26 - f30j | are scheduled to probe the same region of parameter space. The limits of MAMI and 
APEX are those as given in their publications |26l . [27j . Our prediction for the exclusion limit expected in the MAMI 
2012 experiment discussed in section [III Bl is depicted by the dashed curves. The prediction for MESA obtained in 
section llVl is indicated by the dotted (dashed-dotted) curves for the setups with a central scattering angle of 10° (20°). 
Our calculation shows, that the 2012 experiment is well suited to exclude a large region of the parameter space and 
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in particular most of the so-called [g — 2) M welcome band, in which the discrepancy between the experimental and 
theoretical value of the anomalous magnetic moment of the muon (g — 2) M could be due to 7' contribution. 
We propose an experiment for the MESA accelerator under construction at Mainz. The investigated kinematic settings 
will allow for the exclusion of the remaining part of the (<? — 2)^ welcome band that is not probed so far. 



V. CONCLUSIONS AND OUTLOOK 



In this work we have calculated the cross section which are crucial to describe the existing and planned fixed-target 
7' search experiments. A comparison of our calculations with a sample of data taken at MAMI has been performed. 
After applying the leading order QED radiative corrections for the corresponding elastic electron-hadron scattering 
process we find, that our calculations and the data sample are in good agreement. In addition, a calculation of the 
separated spacelike and timelike virtual photon exchange cross sections, each for the direct and exchange term, has 
been performed. This allows us to study the dependence of the background cross section on these contributions. 
Furthermore we find, that it is necessary to include the exchange term into the cross section in order to reconcile the 
data. The exchange contribution is contributing to the irreducible background. 

Using the cross sections obtained in our analysis, we are able to provide predictions for the expected exclusion limits 
for MAMI and MESA. Following our predictions, the experiments at MAMI and MESA will be able to probe the 
entire (g — 2) M welcome band and in addition, increase the existing limits by one order of magnitude. 
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Appendix A: Detailed Cross Section Calculations for 7' Search experiments 

Starting from Eq. one finds by inserting 1 = j d 4 q' 5^ (q' — 1+ — ?_) 

1 d 3 k' d 3 p' d 3 q' 



da = 



4y/(k-p) 2 -m?M 2 (2tt) 3 2^ (2n) 3 2E' p (2tt) 3 2 g'° 
(2n) A S^(k+p-k'-p'-q')- 



{2n) 3 2E^ (2tt) 3 2£; 4 

q'° dq'° 



2tt 



(2TT) i 5^(q'-l + -l-)\M\< 



The (5-functions constrain the three-momenta 

q = /_ + and p' = k — k' — q' , 

which leads to 

1 1 |?'| 2 l/7 I 2 !/" I 2 

da = ^ s J J d\k'\ dn e , d\l+ I dVl + d\l- I dVt_ dq' 2 

128 \k I M (2tt) E p ,E kl E A ,Ei + Ei_ 

S (E + M - E e , - E p , - q'°) S (q'° - E + - E_) ~\M? . (Al) 

V v ' " v ' 

—:6i ='-&2 

The remaining two delta functions can be used to express the energies associated with k' and Z_, by which integration 
over their three- momentum absolute values is performed. Therefore expressions for \k'\ and \l-\ in terms of the 
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remaining quantities have to be found, using 



q' 2 = 2mf + 2E+E- - 2 1 Z \l + ■ Z_ 

0-0= mf] +E+E- - \C 1 1+ ■ L . (A2) 



2 



= :B 



= :A 



This equation can be rewritten as a quadratic equation for |Z_ | which can be easily solved. After adding (B\l_\ — A) 
on both sides of Eq. (jA2[) . squaring the result and using E 2 ^ = |i_ | 2 + mf one finds the two solutions 



AB (AE+r + {E +mi BY - (El mtf 

± \ 7^2 ^2Y?. ■ ( A6 ) 



1 B 2 -E\ y {B 2 -E 2 + ) 2 

The determination of the physical solution can be done by considering the particles as massless. Now the calculation 
simplifies to 

q 12 =l\ + l 2 _ +2|C||Z+|(l-f+-L) 
=0 

<*\C\ = q '\ . . (A4) 

2|Z + |(1 -/+•/_) 

Comparing Eqs. (|A3|1 and (|A4|) . one finds that the solution with "+" corresponds to the physical allowed case. Thus 
it is 



, r , AB j (AE + ) 2 + (E +mi B) 2 - (E 2 +mi ) 2 

{l - 1 ~ & - E% + V (B 2 -E 2 + ) 2 ' (A5j 

with A = -<?' 2 /2 + to 2 and B = 1+ ■ L. 

The calculation of \k'\ is done in a similar way. Since it is not necessary that the four- vectors l + and /_ appear 
explicitely in the following, instead their sum q' 2 = (l + + L) 2 is used where \l- | is symbolic for the result of Eq. (|A5I) . 
Again starting from four-momentum conservation one finds 

& 0= (p + k~ q'f + m 2 - M 2 - 2(E + M - q'°) E e , + 2(k - q') ■ k' \k'\. 

= :D =:F =:G 

An analogous calculation as for |Z_ | then leads to 



Thus one has 



p DG / (mFG) 2 + (DF) 2 - (mF 2 ) 2 

l fc l = -^3^ + \/ (G 2 - F 2 ) 2 ' (A6) 



(E + M — E e , - E v , - q" 



d\k'\ d\k'\ 

\k'\ \k'\-k'-(k-q') 



and 



Ek' E p i 



962 9 tf-E + -EJ) 



(A7) 
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E- q 



/o 



(A8) 



1G 



For the experiments performed at MAMI the detector quantities are given in table [I] The horizontal and vertical 
acceptances are given in a Cartesian reference frame. It is convenient to calculate the cross section directly in the lab 
frame. The lab frame three-momenta of the detected particles depending on these quantities are parametrized by 



\l± 



\j\ + tan 2 56 + tan 2 



' tan 6(p 

\ COS 0o 



cos 4>o + sin 6 
tan 59 
— tan 5(f) sin 4>q 



where 4>o is the central horizontal angle of the detector, 56 is the deviation from the horizontal scattering angle and 
59 is the deviation from the vertical out-of-plane angle. Note that the vertical central angle of the detectors is 0°. 
Integrating over the angles 5<f) and 56 within the limits of the experimental acceptances then leads to the cross section 
Ac. To account for this geometry the cross section has to be multiplied by a Jacobian 



J{56, 56) 



cos 2 56 cos 2 59 (l + tan 2 59 + tan 2 1 



.3/2 



The cross section then reads as 

da 



d\l + \dn+dn-dfl e , dq' 2 
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